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Abstract 

Operating high-intensity discharge lamps in the high frequency range 
(20 —300 kHz) provides energy-saving and cost reduction potentials. How¬ 
ever, commercially available lamp drivers do not make use of this operat¬ 
ing strategy because light intensity fluctuations and even lamp destruction 
are possible. The reason for the fluctuating discharge arc are acoustic res¬ 
onances in this frequency range that are excited in the arc tube. The 
acoustic resonances in turn generate a fluid flow that is caused by the 
acoustic streaming effect. Here, we present a 3D multiphysics model to 
determine the influence of acoustic streaming on the temperature field in 
the vicinity of an acoustic eigenfrequency. In that case a transition from 
stable to instable behavior occurs. The model is able to predict when 
light flicker can be expected. The results are in very good accordance 
with accompanying experiments. 


1 Introduction 

Approximately one fifth of the world production of electric power is consumed 
by 30 to 33 billion electric light sources m- High intensity discharge (HID) 
lamps represent a main part of the artificial light sources. They are used for 
street and shop lighting as well as automobile headlights and other applications. 
A superior color rendering index combined with a sun-like luminance character¬ 
izes HID lamps [4]. These advantages over light emitting diodes will guarantee 
a market for HID lamps in the future. The structure of the low-wattage HID 
lamp, that is investigated in this article, is shown in figure The distance 
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Figure 1: Setup of the investigated HID lamp (Philips 35 W 930 Elite) in hori¬ 
zontal operation position. 


between the tungsten electrodes is 4.8 mm. The tube from polycrystalline alu¬ 
mina (PCA) contains a mixture of argon, mercury and metal halides. A voltage 
is applied to the electrodes to establish a plasma arc of high temperature in¬ 
side the arc tube. In order to minimize material costs of the electronic driver, 
the lamp should be operated in the high frequency range of 20 kHz to 300 kHz 
[5]. However, the periodic ohmic heating excites acoustic resonances in this fre¬ 
quency range. The resonances can induce low frequency plasma arc fluctuations 
of approximately 10 Hz [6] that are visible as light flicker. The acoustic stream¬ 
ing (AS) phenomenon is the link between the high frequency acoustic resonance 
and the low frequency flicker [ 7 ]. 

To investigate the light flicker phenomenon, a three-dimensional stationary 
finite element model of this lamp was developed. Initially, it calculates the 
electric potential, the buoyancy driven velocity field and the temperature field. 
On the basis of the temperature field the acoustic eigenmodes are determined. 
Thereafter, the acoustic pressure field and the AS force are determined for the 
second mode. With the aid of this model, a symmetry breaking transition for 
the flow field was detected [8]. Here, we present results of the backcoupling of 
the superimposed AS and buoyancy force on the velocity and the temperature 
field. 


2 Model 

The finite element model for the numerical investigation of the horizontally 
operated HID lamp comprises the arc tube wall, its filling and both electrodes, 
and it solves several coupled partial differential equations. All simulations were 
executed with COMSOL Multiphysics and MATLAB. 

2.1 Geometry and mesh 

The model of the HID lamp was created with slight simplifications. The coils at 
the electrodes were simplified to cylinders to reduce the number of finite element 
nodes and, therefore, the computing time. The lamp is tilted by 5° against the 
horizontal axis to prevent stability problems due to instable paths related to 


2 













Figure 2: Finite element mesh of the half model. The origin of the coordinate 
system is located in the center between the electrodes. 


the symmetry breaking transition that was observed in a prior investigation 
[8]. The model has to be designed in three dimensions in order to consider 
the upward bending of the plasma arc due to buoyancy. The temperature field 
is symmetrical to the vertical ^- 2 ;-plane (figure [^. Therefore, only one half 
of the physical geometry has to be considered in this simulation step. The 
temperature field is mirrored at the ^-z-plane to map it onto the full geometry 
and to calculate the acoustic modes. 

The finite element mesh (figure used for the simulation of the half model 
consists of approximately 123 000 elements. In the boundary layer of the arc 
tube prism elements were used, whereas the rest of the model is composed 
of tetrahedral elements. The accuracy of the solution was tested with different 
mesh resolutions. The results obtained with the selected mesh differ from results 
with a mesh with 164 000 elements by less than 1 %. The computing time needed 
with the mesh displayed in figure could be reduced by a factor of 10 compared 
to the simulation with the very fine mesh. 

2.2 Temperature Field 

The coupled differential equations of mass, momentum, charge and energy con¬ 
servation were used to obtain the temperature distribution. The differential 
equation of charge conservation 


V • (aS) = 0 (1) 

enables the calculation of the electric field E. The electric conductivity a and 
other material parameters are listed in a previous publication [8]. Because no 
essential differences are expected, a direct current (dc) was implemented instead 
of an alternating current (ac)[^ The inward current density on the electrodes 

^At ac condition the direction of the carrier wave periodically changes. On time average 
the system is symmetric. The dc boundary conditions have been set up accordingly. 
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was obtained by the electrode tip area and the electric current /lamp- With the 
nominal power of the lamp Piamp and the voltage drop between the electrodes 
U the electric current /lamp = Piamp/P was computed. The tube’s wall and the 
symmetry plane were modeled as electric insulators. 

For Newtonian fluids the stationary velocity field u is described by the 
Navier-Stokes equation 


(„--v) 


u = V ' 


— PI + T] ( Wu + fvfl) ] 


(v-s): 


+ / ( 2 ) 


and mass conservation 

V • (pu) = 0 (3) 

(I identity matrix, p density, P pressure inside the arc tube, p dynamic viscos¬ 
ity). Without AS the velocity field is driven by the gravitational force only. Its 
force density is described by 

fi = -Si3pg ( 4 ) 

{ 6 i 3 Kronecker symbol; index 3 corresponds to the Z-direction, g = 9.81 m/s^). 
The boundary condition u = 0 was used at the tube’s wall and the electrodes. 
At the symmetry plane n-u = 0 {n normal vector) was applied and in the center 
the static pressure P = 2.81 MPa was specifiecQ 
The Elenbaas-Heller equation 

V ■ (^-KWTj + pCpU-VT = a\E\^ - Qrad ( 5 ) 

{k thermal conductivity, Cp specific heat capacity, q^ad density of power radia¬ 
tion) determines the temperature field T inside the arc tube, of the wall and of 
the electrodes. On the external surface of the wall heat is emitted by radiation. 
This is described by the boundary condition n • (—/^VT) = ease (Pamb ~ 

(e emissivity, asB Stefan-Boltzmann constant, Tamb ambient temperature). On 
the symmetry plane the emissivity is set to zero (no radiation). At the elec¬ 
trode grounds a constant temperature of T = 1450 K was defined. The acoustic 
pressure field in the arc tube can be calculated with the aid of the resulting 
temperature field and the power density of heat generation. 


2.3 Acoustic Streaming 

For the calculation of the AS force the acoustic pressure in the arc tube has to 
be calculated first. To obtain the acoustic pressure p, the Helmholtz equation 

^ / 1 ^ \ (jj‘^ 'y — 1 

V • -Vp H- yp = iuj H ( 6 ) 

\p J pc^ pc^ 

has to be solved (c speed of sound, 7 heat capacity ratio). The source term of 
the Elenbaas-Heller equation P = cr|Pp — q^^d describes the power density of 

^To preserve model symmetry, the center point was chosen. The position of the static 
pressure point is of no importance. 
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heat generation. The arc tube filling is treated as an ideal gas with the speed 
of sound field 


c = 



( 7 ) 


The wall of the arc tube and the electrodes are treated as sound hard. The 
solution of the Helmholtz equation is expressed by an acoustic mode. It is 
described elsewhere, how the pressure amplitude, that takes the loss mechanisms 
of heat conduction and viscosity on the surface as well as in the volume into 
account, is calculated isHn]. 

The AS force arises when a standing sound wave in the arc tube of the HID 
lamp is excited [7]. In addition to an oscillatory motion, an averaged fluid mo¬ 
tion is induced. The motion is a consequence of the viscous force that develops 
due to the shear viscosity of the fluid particles. It generates vortex-like motions 
that are distinguished by their location: Outer (or Rayleigh) streaming occurs 
outside the boundary layer and inner streaming occurs inside the boundary layer 

m- 

The AS force is implemented into the source term of the Navier-Stokes equa¬ 
tion. Instead of the gravitational force only, the force density is now 


fi 


dpVkVi 

dxk 


^ispg 


( 8 ) 


{v sound particle velocity). In this equation Einstein’s summation convention 
is used and the overbar indicates time averaging. In the case of time harmonic 
waves the force density simplifies to 


1 dpvkVi 

2 dxk 


^i3pg- 


(9) 


V denotes the amplitude of the sound particle velocity. It can be obtained from 
the acoustic pressure [13]: 


V (r, uj) = :—Vp (r, uj ). 
loop 


( 10 ) 


The velocity drop in the boundary layer had to be considered additionally be¬ 
cause the mode calculation does not incorporate this condition. Therefore, the 
sound particle velocity was multiplied with a factor that drops in the boundary 
layer from one to zero nans]. The implementation of this function into the 
model has previously been described [8]. 


2.4 Simulation procedure 

The calculation of the temperature, the acoustic pressure and the AS force 
field have to be coupled in order to get the stationary temperature field at a 
certain frequency. Since the discharge arc is excited at this frequency, it is called 
excitation frequency. 
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Initially, the velocity and temperature fields as well as the acoustic mode 
and eigenfrequencies are calculated with the gravitational force only. On the 
basis of the acoustic modes the pressure amplitude is calculated. With the aid of 
the amplitude, the AS force at a specified excitation frequency is computed. In 
order to include the back coupling of the AS field on the other fields, a recursion 
loop was implemented. From the second iteration on, the fluid flow is driven by 
buoyancy and the AS force. This leads to a different velocity and temperature 
field compared to the first step. The acoustic eigenmode and the resonance 
frequency change and lead to an altered AS force that serves as input for the 
next step. This recursion procedure is continued until the results of the current 
and the previous step coincide within a certain limit. The resonance frequency 
is used as a convergence criterion. Once the resonance frequency changes by 
less than 5 Hz, the solution is considered converged. 

The fields obtained in this way are used as initial condition for the next 
excitation frequency. Since we are mainly interested in results near a resonance, 
the excitation frequency is increased when the resonance frequency is higher and 
decreased when the resonance frequency is lower than the excitation frequency. 
The step size is adopted to the difference of both frequencies. The closer the 
resonance and the excitation frequency are, the higher the pressure amplitude 
and AS force. Because the AS force depends quadratically on the acoustic parti¬ 
cle velocity (see Equation]^, the force increases significantly near the resonance 
frequency. Hence, a very small step size of the excitation frequency has to be 
chosen. 

3 Experimental setup 

In the experiment the operating voltage of the lamp consists of two components. 
The square-wave part with a carrier frequency fc of 400 Hz is used to operate 
the lamp at stable conditions. Flicker of the discharge arc is excited by a high 
frequency sinusoidal part. This excitation frequency fe is varied in the range 
of 41kHz to 44kHz. The combined voltage input is described by (cje = 27r/e, 
ujc = 27r/c) 



2 n < ujct < 2n + 1 
2 n + 1 < uJct < 2n + 2 


n e No. (11) 


The amplitude A is adjusted to the nominal power of the lamp. The ratio of 
the square-wave to the sinusoidal-wave voltage amplitude is described by the 
modulation depth a. In the experiments it is set to a constant value that is 
just high enough to stimulate arc flicker. A more detailed description of the 
experimental setup can be found elsewhere [16]. 

Before the measurements the lamp is operated at stable conditions {a = 0) 
for at least 20 minutes to reach a stationary state. Then, the modulation depth 
is set to a constant value and, simultaneously, the excitation frequency is set 
to its initial frequency some kHz off the resonance frequency. This initial value 
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is either 41kHz when the frequency is afterwards increased or 44 kHz when 
the frequency is decreased. These operating conditions are kept constant for 
20 s. During this time the voltage drop between the electrodes is measured 40 
times. At the initial step no acoustic resonance is excited so that the mean 
value of the measured voltages can be defined as the basic voltage. In the 
next step the excitation frequency is shifted by 50 Hz into the direction of the 
resonance frequency. To prevent lamp failures, which are caused by too high arc 
tube temperatures or by temperature oscillations of the arc tube, the voltage 
measurements are from then on used to control the experiment. If one of the 
40 measurement exceeds the basic voltage by 8 V or if the voltage fluctuates by 
more than 1.5 V, the experiment is terminated. If these criteria are not fulfilled, 
the experiment is continued until the excitation frequency reaches 41 kHz and 
44 kHz, respectively. 


4 Results 

4.1 Acoustic modes and eigenfrequencies 

Initially, the 5° tilted model without AS was solved. The simulation computes a 
lowest acoustic eigenfrequency of 33 100 Hz which is not apparent in the exper¬ 
iment. At 47 750 Hz and 48130 Hz two neighboring azimuthal modes appeal]^ 
In this article the simulations, which include the backcoupling of AS on all other 
fields, focus on the mode at 47 750 Hz. The experiment shows the influence of 
both eigenfrequencies. 

The excitation frequency was set to 47 750 Hz and the recursion started. 
Once convergence was reached, the eigenfrequency shifts to 47400 Hz due to 
AS. As described above, the excitation frequency for the next step has to be 
decreased by 50 Hz. The resulting eigenfrequency is 40 Hz below the previous 
value at 47360 Hz. Reducing the excitation frequency again by 50 Hz further 
decreases the difference between excitation and eigenfrequency so that the sim¬ 
ulation of the next step is conducted closer to the resonance. This results in a 
stronger AS force. Continuing in this way the simulations were successful down 
to an excitation frequency of 46 313 Hz where the difference to the resonance 
frequency diminishes to 55 Hz. At even lower frequencies the simulation did not 
converge. We assume that no time-independent solution exists at these frequen¬ 
cies so that stationary simulations are inappropriat^ To test this assumption, 
we plan to perform a linear stability analysis in the near future. 

4.2 Voltage drop 

A quantity, which is easy to determine in the simulation as well as in the ex¬ 
periment, is the voltage drop between the electrodes. In order to demonstrate 

gravity is switched off, the two modes are degenerate. 

similar situation is known from the Taylor-Couette system, where the stationary Taylor 
cells undergo a Hopf bifurcation to the time-periodic wavy vortex state El. 
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Figure 3: Simulation result of the voltage drop between the electrodes over the 
excitation frequency. The blue crosses mark voltage drops at decreasing exci¬ 
tation frequency. The red circles indicate voltage drops at increasing excitation 
frequencies. 


the reliability of the FE model, the voltage of the simulations and experiments 
are compared. 

Figure shows the simulation results of the voltage as function of the ex¬ 
citation frequency. The two curves correspond to the behavior of the voltage 
for increasing and decreasing excitation frequency, respectively. When the fre¬ 
quency is decreased from 47.8 kHz to 46.3 kHz, a significant voltage increase 
of 7V occurs. A further lowering of the frequency leads to a voltage jump to 
values coinciding with the values of the curve representing the increasing fre¬ 
quency simulation. At increasing excitation frequency the voltage stays constant 
at about 84 V for the frequency range from 45.0 kHz to 47.0 kHz. Then the volt¬ 
age increases by about 1V, before it jumps to the upper curve. The data clearly 
demonstrate the appearance of a hysteresis, which is expected in a nonlinear 
system. The situation is similar to the one of a Duffing oscillator: The reso¬ 
nance frequency depends on the amplitude. If the oscillation is strong enough, 
the equation, which couples the amplitude and the excitation frequency, is cubic 
and its solution includes three real roots in a certain frequency range. The peak 
of the response curve either tilts to low or high frequencies and overlaps the 
response curve describing the amplitude far off the resonance [18]. 

Figure]^ shows the results of the experimental voltage measurements over the 
excitation frequency. The frequency deviation between the experimental and the 
simulation result of 4 kHz is related to the simplified geometry and uncertainties 
in the material properties [Ellin]. Additionally, some boundary conditions, e.g. 
the constant temperature at the electrode ground or the uniformly distributed 
electric current density on the electrodes, only approximate the conditions in 
the lamp. For a decreasing excitation frequency the voltage increases up to 96 V 
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Figure 4: Experimental results of the voltage drop between the electrodes over 
the excitation frequency at a modulation depth of 2%. Red circles indicate 
measurements at increasing and blue crosses at decreasing excitation frequency. 


at 42.4 kHz. Here, the experimental results show a voltage jump like the one 
observed in the simulation. This jump has also been observed in high pressure 
sodium lamps [20]. The significant voltage increase is related to the excitation 
of an acoustic mode. In contrast to the simulation, in which only one of the 
neighboring acoustic modes is considered, figure shows repercussions of both 
modes. The small voltage peak at the frequency of 43.8 kHz indicates the second 
mode. 

The simplified geometry and uncertainties in the material properties are also 
responsible for the different voltages: In the simulation the voltage far off an 
acoustic resonance is 83.8 V and in the experiment approximately 91.5 V. Apart 
from these differences in the absolute values the graphs in figure and figure 
show excellent agreement. 

The voltage jump at an increasing excitation frequency is 2.5 V compared to 

2.3 V in the experiment. At a decreasing excitation frequency the voltage in the 
experiment increases by 5.1V and in the simulation by 8.3 V. The frequency 
range of the hysteresis (voltage peak at an increasing to the voltage peak at 
a decreasing excitation frequency) is 840 Hz in the simulation and 500 Hz in 
the experiment. It would be interesting to measure the hysteresis at different 
modulation depths to determine its influence on the voltage jump and the jump 
frequencies. 

4.3 Temperature and velocity field 

Figure [fright shows the temperature field resulting from the simulation at an 
excitation frequency of 46 313 Hz, the frequency, where the voltage jump in the 
blue curve of figure [^appears (decreasing frequency). A major impact of AS on 
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Figure 5: Temperature distribution in the ^-z-plane at different excitation fre¬ 
quencies while decreasing the excitation frequency. Left: Far off resonance 
(47 750 Hz). Right: Near resonance (46 313 Hz). The black line indicates the 
location of maximal temperature as a function of the ^-coordinate. 


the temperature field can be observed, when comparing it to the temperature 
field obtained at 47 750Hz (figure left). This frequency is far off an acoustic 
resonance and the AS force accordingly small. 

While approaching the resonance, the maximum plasma temperature, lo¬ 
cated at the hot-spots in front of the electrodes, only increases slightly. Simul¬ 
taneously, the temperature in the lower part of the arc tube, especially in the 
vertical center plane of the arc tube, considerably increases. These higher tem¬ 
peratures are caused by a flow that is directed in negative ^-direction. The flow 
in the upper part is directed into the positive 2 :-direction and leads to a steeper 
temperature gradient near the upper wall. 

To compare the temperature field of the simulation with experimental re¬ 
sults, the light intensity distribution was recorded by a camera [16]. The left 
part of figure shows the intensity distribution at stable operation conditions. 
The buoyancy force bends the plasma arc upward and AS is negligible. The 
picture resembles the left part of figure where the temperature distribution 
far off the resonance is depicted. 

After tuning the parameters to unstable lamp operation, AS emerges. The 
intensity distribution on the right shows the discharge arc when it moves towards 
the lower part of the arc tube. Immediately afterwards, the discharge arc started 
to rotate around the ^-axis and finally the lamp exploded. It highlights that AS 
acts against buoyancy. The pictures strongly resemble the right part of figure 
from the simulation. 

Figure shows the simulated flow fields corresponding to the temperature 
fields depicted in figure At an excitation frequency far off the resonance the 
fluid flow essentially represents the purely buoyancy-driven flow field. The flow 
between the electrodes is directed upwards and has its maximum in this region. 
Instead of a maximal velocity of 0.09 m/s in case of the purely buoyancy-driven 
flow, the maximal velocity is 0.13 m/s. The flow field is not fully symmetric 
with respect to the ^-z-plane due to the tilting of the lamp. 

Near the resonance the AS dominates the buoyancy flow and the asymmetry 
caused by the lamp tilting is not visible any more because only the buoyancy 
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Figure 6: Measurements of the light intensity. The left part shows the plasma arc 
at stable operation conditions. The right part shows the plasma arc immediately 
before it starts to rotate violently around the lamp axis. The black line indicates 
the location of maximal intensity as a function of the ^-coordinate. 



Figure 7: Fluid flow in the ^-z-plane at different excitation frequencies while 
decreasing the excitation frequency: Left: Far off resonance (47 750 Hz). Right: 
Near resonance (46 313 Hz). 
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force is affected by the tilting of the lamp. In the upper part of the arc tube 
two additional vortices develop (right: clockwise; left: counterclockwise). In 
the lower part the flow points in the direction of the negative 2 :-axis. The flow 
continues in the x-direction, which is not visible in the figure. The maximal 
velocity is 0.42 m/s whereas in the simulation without backcoupling it is 0.61 m/s 
[8]. The reason for the difference is that in case of the coupled simulation the 
lamp is excited 55 Hz off the resonance. The AS force is too weak to cause the 
symmetry breaking transition, which was observed in the uncoupled model [8]. 

5 Conclusions 

In this work a stationary 3D finite element model is used to simulate the back- 
coupling of acoustic streaming on the temperature field inside a high-intensity 
discharge lamp. At the investigated acoustic eigenfrequency a hysteresis was 
detected in the simulation and confirmed by accompanying experimental inves¬ 
tigations. The voltage drop when decreasing the excitation frequency shows a 
different behavior than in case of an increasing excitation frequency. The sim¬ 
ulation results show a severe impact of the acoustic streaming on the voltage 
drop, the temperature field and the flow field. In the simulation as well as 
in the experiments the voltage drop at the investigated resonance significantly 
increases compared to a voltage drop far off the acoustic resonance. 

AS causes the maximal velocity in the arc tube to become more than 4 times 
higher compared to the purely buoyancy-driven flow. A buoyancy dominated 
flow leads to an upward bending of the discharge arc. Acoustic streaming on 
the contrary induces a straight arc between the electrodes. These results are 
consistent with measurements of the light intensity distribution. The results 
will help to understand the physical processes, that lead to flickering of the 
discharge arc inside high-intensity discharge lamps. They enable development 
of new lamp systems with improved energy-efficiency. 
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